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Preface 


The random fluctuations generated by the turbulent boundary layer are a dominant source 
of interior cabin noise at subsonic as well as supersonic speeds; added treatment to man- 
age this noise usually is a significant weight consideration. Fluctuating surface pressure 
measurements were conducted on flight tests of the Tupolev 144LL Flying Laboratory in 
late 1997 and early 1998 in Russia. These flight tests used pressure sensors mounted in 
cabin window blanks in different locations down the length of the fuselage. The recorded 
data were returned to the U.S. for analysis by NASA and Boeing. 

The literature on turbulent boundary layer pressure fluctuations provides several empirical 
models which were compared to the measured TU-144 data. The model developed by 
Efimtsov showed the best agreement. Adjustments were made to improve its agreement 
further, consisting of the addition of a broad band peak in the mid frequencies, and a minor 
modification to the high frequency roll-off. The adjusted Efimtsov predicted and measured 
results are compared for both subsonic and supersonic flight conditions. Measurements in 
the forward and middle portions of the fuselage have better agreement with the model than 
those from the aft portion. For High Speed Civil Transport supersonic cruise, interior levels 
predicted by use of this model are expected to increase by 1 to 3 dB due to the adjust- 
ments to the Efimtsov model. 

The space-time cross-correlations and cross-spectra of the fluctuating surface pressure 
were also investigated. This analysis is an important ingredient in structural acoustic mod- 
els of aircraft interior noise. Once again the measured data were compared to the pre- 
dicted levels from the Efimtsov model. No adjustments were made to the cross-correlation 
portion of the prediction model. 

The auto-spectral model presented in this report can be used as improved input data to 
structural acoustic prediction models. However, further investigation of the cross- 
functional model is recommended to fully utilize the information available from the TU-144 
measured data. 
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1. Introduction 


The motivation for this work stems from the need for an aerodynamic loads model for predicting 
passenger aircraft interior noise, at subsonic as well as supersonic cruising speeds. The turbulent 
boundary layer on the fuselage causes structural vibrations which radiate random noise into cock- 
pit and passenger cabin. Simple interior noise prediction models use auto spectra of the pressure 
fluctuations. More complex models require space-time cross-correlations (or cross-spectra) of 
pressure fluctuations. 

This work was performed at the Boeing Company under NASA Contract 20220, High Speed Re- 
search Program, Phase II, Task 27, Structural Acoustics. The objective was the validation of fluc- 
tuating pressure turbulent boundary layer models or model, both for single signals (auto spectra) 
and signal pairs (cross spectra). 

The data set underlying this study came from measurements during flight tests of the Tupolev 144 
Supersonic Flying Laboratory in late 1997 and early 1998 [Ref. 10], 

2. Models 

Several models available in the literature are discussed and then used to compare the present 
data with their predictions. A model by Smol’yakov-Tkachenko [Ref. 17] is not used because it 
was built with marine applications in mind and is not well suited to high speed flight [Ref. 11], As 
discussed in the Sections below, several authors have built different models for power spectra 
(auto spectra) at a single point of the surface. In those references, models for space-time cross- 
correlation or frequency-wave number cross spectra are either not discussed, or are cast in only 
one form (exponential decay with spatial separation) which is discussed in Section 2.5 (Efimtsov). 
The equations presented here are meant to convey the flavor and complexity of the models. If the 
reader wishes to actually implement any of these, please refer to the complete references. 

2.1 Robertson [Ref. 1] 

This prediction of the turbulent boundary layer fluctuating pressure power spectrum is derived from 
an empirical evaluation of sub-sonic and super-sonic wind tunnel experimental results. The normal- 
ized equation for the frequency spectrum is a function of an estimate of the mean squared pres- 
sure and the Strouhal number. 
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2.2 ESDU [Ref. 2] 

The predicted power spectrum is based on an analysis of flight test results. Fluctuating pressure 
levels, flow field, and boundary layer parameters from flight testing is curve fitted to determine the 
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predicted octave and third octave band levels. The prediction is a function of: Mach number, 
Strouhal number, and Reynolds number. The curve fitting equations appear in Section 10 (Appen- 
dix A). 


2.3 Laganelli [Ref. 3] 

Predicted power spectra are derived from experimental and fluid dynamic principles. Experimental 
data includes: AFWAL Mach 3 facility, Blake, Laganelli, Raman, Coe, and Bull (these do not ap- 
pear to include flight data). The equation is an expansion of Robertson's equation, but adds fric- 
tion coefficient considerations: 
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where F c is a transformation from compressible to incompressible flow states: 
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2.4 Goodwin [Ref. 4] 

This fluctuating pressure predicted spectrum was derived from previous prediction equations and 
flight test data on three super-sonic test aircraft: XB-70, A3J, and Concorde. A very small portion 
of the XB-70 data and even less of the A3J data was usable for this type of prediction due to a 
number of experimental errors. Additionally, the Concorde data, which used Kulites in window 
blanks very similar to those used for the Tu144 testing, is somewhat suspect due to the experi- 
menter not paying close attention to the sensitivity of the Kulites to the flushness of the sensor's 
surface relative to the surrounding surface (discussed further in Section 3.1.3). The Kulite meas- 
urements varied as much as 9 dB between sensors in the same window blank. This test also ex- 
perienced trouble duplicating results for repeated test conditions. Goodwin adjusted the Robertson 
and Laganelli predictions to account for the measured in-flight data using the following equations: 
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2.5 Efimtsov [Ref. 5 and 13] 

Professor B. M. Efimtsov of TsAGI (Central Aerohydrodynamic Institute in Moscow, Russia) up- 
dated the existing model of turbulent boundary layer fluctuating pressure levels. This new model 
utilized data from the low and high speed TsAGI wind tunnels, Tu144, and what appears to be a 
TU-22 (twin engined, supersonic military aircraft) in-flight measurements. Details of the flight con- 
ditions and sensor locations are given; they appear to be similar to those measured on the Tu144 
by NASA-Boeing. However, no actual flight data is presented for direct comparison to the NASA- 
Boeing tests. 

The analytical derivation of the TsAGI predicted power spectrum references Laganelli's calcula- 
tions; the principal independent variables are Strouhal number and Reynolds number. Efimtsov's 
equation is: 

® M = 2 noUlp 1 8 — — ~iy , 

(l + 8 a 3 Sh 2 y 3 +cc]3Re( s y Re J 


fi = 



,a = 0.01 , 


Sh = -jj- (Strouhal Number) 


The model for cross-spectral density assumes that dependencies on frequency and spatial separa- 
tion can be mathematically separated as follows: 
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where subscript 1 refers to the longitudinal (flow-wise) direction, and subscript 2 to the lateral 
(cross flow along the surface) direction. ^ and £ 2 are separation distances between points on the 

surface, and A l and A 2 are the correlation length scales which strongly depend on frequency (see 
Ref. 13 for details). 


3. Boundary Layer Data 

The reader is referred to Ref. 10 for details on how TU-144 data were acquired. Figure 1 shows the 
flight conditions for which data was acquired. Due to lack of time for this study, only a representa- 
tive subset of the data was utilized. 
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Altitude, 


TU-144 FLIGHT ENVELOPE 



Indicated Air Speed, km/hr 

Figure 1 TU-144 Flight Envelope and Experiment 2.1 Test Points 


Figure 2 shows the layout and nomenclature of the transducers in the window blanks and in the 
fuselage skin. 
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N1.3 N1.1 N1.5 
N1.4 


Figure 2 Location and Identification of Window Blanks (WB) Instrumented with Kulite Pressure 

Transducers (N1.1 through N7.5, SI, S2) 
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3.1 Data Analysis 


3. 1. 1 Data Reduction 

Pressure level spectra were calculated from one minute samples recorded for each transducer and 
for each flight condition. These spectra were calculated both on a narrow band and a third octave 
band basis. The narrow bands have a uniform constant bandwidth of 2.44 Hz; spectral levels re- 
ported here were corrected to a 1Hz bandwidth. The third octave band levels were also corrected 
to a 1 Hz bandwidth for direct comparison between the narrow band and the third OB and to com- 
pare levels to some of the published results from previous investigations. 

All signal analysis was carried out using the MATLAB software package on an HP-UX work station. 
Analysis scripts were provided by Dr. Stephen A. Rizzi of NASA Langley Research Center. 

3.1.2 Data Normalization 

Two forms of data normalization were investigated to help compare between channels and be- 
tween flight conditions, and for comparing data to previous tests. The first normalization is what is 
seen in most of the literature on turbulent boundary layer spectral shapes [Ref. 6 and 7] 

0(co)U 0 cod* 

— VA; — VS. 

q~8 U 0 

where U 0 is the free stream velocity, q 0 is the dynamic pressure, and 8* is the displacement 
boundary layer thickness. Efimtsov [Ref. 5] suggests a slightly different normalization for a better 
collapse of the TsAGI measured data: 

0)8 

, / vs. 

p 2 U^8 U 

r w T T 

where U x is the friction velocity defined by the quotient of wall shear stress over fluid density at the 
wall p w , and 8 is the (velocity deficit) boundary layer thickness. We usually use the Efimtsov nor- 
malization in our work. 

3.1.3 Kulite Flushness 

An exhaustive wind tunnel test was performed by Professor Efimtsov at TsAGI [Ref. 8] to deter- 
mine the effect sensor flushness would have on the measured results at both subsonic and super- 
sonic speeds. For measured flushness ranging from recessed to protruding, deviation from ex- 
pected third octave band levels were reported. Variations of less than 2 dB can be expected for 
sensors deviating less than 20x1 0" 6 m from flush, the flushness deviation being measured in the 
stream-wise direction. The Tu144 flight test Kulite flushness measurements are discussed in Ref. 
10. Only four of the Kulites measured this much of a deviation from flushness. However, when 
these four Kulites were compared to Kulites in the same window blanks (Section 3.1.4), only one of 
them (N4.1) varied significantly from the average measured SPL within that blank; data from Kulite 
N4.1 was not used in the auto spectrum analysis in this study. 

3 . 1.4 Outlying Fluctuating Pressure Levels 

Within window blank comparisons were conducted by plotting the third octave band levels for all of 
the sensors within each window blank for 22 selected steady-state flight conditions. Channels that 
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showed significant deviation from the other channels within the same window blank were noted as 
deviant and were not used for further auto-spectral analysis (Table 1). Causes for the deviation in 
measured SPL may included: line noise, damaged sensors, sensor flushness, localized flow dis- 
turbances. 


Table 1 Outlying Fluctuating Pressure Measurements. 


Flight & 
Run # 

Altitude 

(kft) 

Mach # 

Window 
Blank 1 
Only 

Window 
Blank 4 
Only 

Window 
Blank 7 
Only 

Central 
Transducer 
on all 7 
Window 
Blanks 

Bad Channel #'s 

f9r17 

16.4 

0.648 

OK 

4.1, 4.9 

OK 

OK 

5 

26 



f9r21 

16.4 

0.877 

OK 

4.1, 4.9 

OK 

7.1 

5 

16 

26 


f9r23 

13.1 

0.563 

OK 

4.1, 4.9 

OK 

OK 

5 

26 



f10r29 

45 

1.6 

OK 

4.1, 4.4, 
4.5 

7.5 

4.1 

5 

8 

9 

30 

f10r31 

46.6 

1.6 

OK 

4.1, 4.5 

7.5 

4.1 

5 

9 

30 


f10r33 

55.4 

1.95 

OK 

4.1 

7.5 

4.1 

5 

30 



f10r34 

55.4 

1.95 


4.1, 4.3 


4.1 

5 

7 



f10r35 

56.8 

1.95 

OK 

4.1, 4.3 

7.5 

4.1 

5 

7 

30 


f10r36 

56.8 

1.95 


4.1, 4.3 


4.1 

5 

7 



f10r37 

56.8 

1.95 


4.1, 4.3 


4.1 

5 

7 



f 1 1 r40 

52.5 

1.8 

OK 

4.1, 4.5 

OK 

4.1 

5 

9 



f 1 1 r42 

54.5 

1.7 

OK 

4.1, 4.3, 
4.5 

OK 

4.1 


H 

9 


f 1 1 r44 

50.9 

1.7 

OK 

4.1, 4.5 

OK 

4.1 

5 

9 



f 1 1 r48 

41.3 

1.4 

OK 

4.1, 4.5 

OK 

4.1 

5 

9 



f 1 1 r50 

41 

1.2 

1.2 

4.1 

OK 

4.1 

5 

22 



f 1 1 r53 

16.4 

0.75 

OK 

4.1 

OK 

4.1 

5 




f17r66 

53.1 

2 

OK 

4.3 

OK 

4.1 

5 

7 



f17r68 

55.4 

2.05 

OK 

4.3 

OK 

4.1 

5 

7 



f17r70 

53.8 

1.75 

OK 

4.1, 4.3, 
4.5 

OK 

4.1 


7 

9 


f 1 7r71 

49.2 

1.6 

N/A 

4.3 

N/A 

OK 

m 





Single sensors located on window blanks by themselves (blanks 2,3,5, and 6) and the two sensors 
on the skin (SI and S2) were more difficult to compare for outlying SPL. Normalized levels plotted 
for sensors down the length of the plane could be used to some extent. However, the uncertainties 
in the data were too great for meaningful interpolations so that data from single sensors were not 
included in the auto-spectral analyses in this report. 
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4. Data Correction 


A correction was added to all auto-spectra of boundary layer surface pressure fluctuations. It ac- 
counts for the lack for spatial turbulence scale resolution of the transducer at high frequencies, ac- 
cording to Corcos [Ref. 14]. Ideally, one would follow the theory presented in Ref. 14 and develop 
the correction formula based on the transducer’s actual geometry as shown in Figure 3. This was 
attempted but abandoned due to mathematical complexity and because an assumption made in 
the theory was not satisfied which says that the sensing areas on the transducer face contribute to 
the total instantaneous measured pressure without phase shift. In the present case, the 12 small 
circles in Figure 3 are just through holes to a small volume behind the face containing a single 
strain gage to which each hole contributes with unknown phase shifts. Instead, a circular trans- 
ducer sensing surface was assumed in order to use the calculated results in Ref. 14. 


Then, the crucial question is what transducer radius to use in the correction procedure. Efimtsov 
[Ref. 8] stated that the transducer he investigated in the wind tunnel behaved with regard to a Cor- 
cos correction as if its diameter was about 8 mm, i.e. double its outer diameter. However, when 
we used a radius of 4mm for the Corcos correction, the high frequency behavior of the corrected 
spectra became unreasonable: instead of a more or less constant roll-off slope on a log-log plot, 
the spectrum started to curve upwards again at higher frequencies. Instead, Philippe Spalart [Ref. 
15] suggested that an equivalent radius be determined on the basis of calculating root-mean- 
square distances from points on the sensing device to the center. With a ring, the root-mean- 
square equals the radius of the ring. With a disk, it is the disk radius divided by square root of two. 
Therefore, we multiplied the radius of the circle on which the sensing holes are located (the “ring”) 


by V2 resulting in 


0.104 




inches = 1.87 


mm . 


This resulted in reasonable corrections. An ex- 


ample is shown in Figure 4. 
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Figure 3 Geometry of Face of Kulite Transducer XCS-190 
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Corcos Correction Comparison, M=1.95, Alt=54000 ft 



5. Comparisons of Measured Data with Models 

5.1 Auto-Spectra 

Figure 5 and Figure 6 (absolute levels) and Figure 7 (normalized) show comparisons between all 
the models and the measured data, for window blank 4 (in the middle of the fuselage). The 
Efimtsov and Goodwin models are the most accurate for predicting the measurements on the 
Tu144. Figure 8 and Figure 9 contain comparisons of TU-144LL data with Efimtsov and Goodwin 
predictions for window blanks 1, 4, and 7. These data are typical of what we saw: at subsonic 
speeds, the Efimtsov model predicts well, and at supersonic speeds both the Efimtsov and Good- 
win models provide good estimates. All models but Efimtsov’s predict a steeper high frequency roll- 
off than the TU-144LL data (which has been corrected for transducer size). Efimtsov's model was 
chosen for further adjustment, as it had the best overall match to the measured data. 
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Predicted and Measured SPL Data, Mach 0.74, Alt=16400 ft 
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90 


Figure 5 Comparison of the Five Investigated Models, Robertson, ESDU, Laganelli, Good- 
win, and Efimtsov, and Measured Data on Tu144 (M 0.74, 16,400 ft, Window Blank 4). 
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dB, (ref: 2x10'°Pa) 


Predicted and Measured SPL Data, Mach 1 .95, Alt=55400 ft 



Figure 6 Comparison of the Five Investigated Models, Robertson, ESDU, Laganelli, Good- 
win, and Efimtsov and Measured Data on Tu144 (M 1.95, 55,400 ft, Window Blank 4). 
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Normalized Predicted and Measured Data, Mach 1 ,95,Alt=55400 ft 



V 

Measured Window Blank 4 


Predicted Efimtsov 

— 

Predicted Goodwin 

— 

Predicted ESDU 


Predicted Laganelli 


Predicted Robertson 


co8/U 

X 

Figure 7 Normalized Comparison of the Five Investigated Models, Robertson, ESDU, La- 
ganelli, Goodwin, and Efimtsov and Measured Data on Tu144 (M 1.95, 55,400 ft, Window 

Blank 4). 
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Predicted and Measured SPL Data, Mach 0.74, Alt=16400 ft 
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Measured Window Blank 7 
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Goodwin Predicted WB 7 


Efimtsov Predicted WB 7 
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Measured Window Blank 4 
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Efimtsov Predicted WB 4 
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Measured Window Blank 1 
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Goodwin Predicted WB 1 


Efimtsov Predicted WB 1 
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Figure 8 Closer Comparison of Efimtsov and Goodwin Models and In-Flight SPL Data (M 

0.74, 16,400 ft, Window Blanks 1, 4, and 7) 
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Predicted and Measured SPL Data, Mach 1.95, Alt=55400 ft 
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Figure 9 Closer Comparison of Efimtsov and Goodwin Models and In-Flight SPL Data 
(M 1.95, 55,400 ft, Window Blanks 1, 4, and 7) 
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5.2 Cross-Functions 

All measured cross-functions presented here were computed using the entire 1 minute time history 
records available in the TU-144LL data. The cross-correlations use the entire frequency band 
available: from about 2 Hz to 11.2 KHz. The cross-spectra were calculated on a narrow band ba- 
sis (2.44 Hz) and then averaged over one third octave bands, corrected to a 1 Hz bandwidth. 

A correction for transducer size was not applied (analogous to the Corcos [Ref. 14] correction for 
auto spectra). This would have required applying digital filters to amplify the high frequency portion 
of the two input time histories being cross-analyzed. 

5. 2. 1 Cross-Correia tions 

Cross-correlations are defined by 

1 T r 

lim^rj p(x,t)p(x+$,t+T)dT 
<p(x, %,r) = r ~ r n r , 

b( x )],J^( x+ ^L 

where p is the fluctuating surface pressure, xis a location on the fuselage, ^ = (<^^ 2 ) is a spatial 

separation vector, t is the time, and r is a time delay. In the following text the location on the fuse- 
lage is omitted and implied by context. 

All the cross-correlations exhibit expected behavior: 

• For longitudinal separations (stream-wise) the peak occurs at a time delay corresponding to the 
convection velocity, and the magnitude at the peak decreases with increasing separation. The 
cross-correlations are well behaved and show discernible correlations even at the largest sepa- 
ration distance of 1741 mm (5.7 feet). The convection velocity can be calculated by dividing 
the separation distance into the time delay of the correlation peak. The result is shown in 
Figure 19. 

• For lateral separations, the peaks occur close to zero time delay. 

• For equal separation distances, the magnitude for longitudinal separations is much greater 
than for lateral separations. 

The angle of local surface stream lines was not known; therefore, no corrections have been ap- 
plied for the flow not being aligned with the transducers. For some lateral separations, the peak 
occurs slightly before zero time delay which may be explained by the flow direction not being per- 
pendicular to the line between the transducers. 

5. 2. 1.1 Cross-Correlations at Front Fuselage 

Figure 10 and Figure 11 show cross-correlations with longitudinal (stream-wise) separation dis- 
tance as a parameter, for transducers on window blank 1. 

Figure 12 and Figure 13 show cross-correlations, with lateral (cross stream) separation distance as 
a parameter, for transducers on window blank 1. 
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Figure 10 
Cross- 
correlations - 
Longitudinal, 
Window 
Blank 1, Wide 
Band, Mach 
0.74 


Figure 11 
Cross- 
correlations - 
Longitudinal, 
Window 
Blank 1, Wide 
Band, Mach 
1.95 
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Figure 12 
Cross- 
correlations - 
Lateral, Window 
Blank 1, Wide 
Band, Mach 
0.74 


0.2 - 




Figure 13 
Cross- 
correlations - 
Lateral, Window 
Blank 1, Wide 
Band, Mach 
1.95 
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5. 2. 1.2 Cross-Correlations at Center Fuselage 

Figure 14 through Figure 18 show cross-correlations which are functions of time delay, with sepa- 
ration distance as a parameter, for transducers on window blanks 3 through 5, and between them. 

The increase in peak value of the correlation going from 60mm to 120mm spacing, as shown in 
Figure 14 and in Figure 16, can be explained by the inaccurate auto spectrum measured on trans- 
ducer N4.1 (discussed in Section 3.1.3). 



Figure 14 
Cross- 
correlations - 
Longitudinal, 
Window Blanks 
3/4/5, Wide 
Band, Mach 
0.74 
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Figure 15 
Cross- 
correlations - 
Longitudinal, 
Window Blanks 
3/4/5, Wide 
Band, Mach 
0.88 


Figure 16 
Cross- 
correlations - 
Longitudinal, 
Window Blanks 
3/4/5, Wide 
Band, Mach 1.2 


20 


«§!■ ty«0.0) <K£ r t)/i|)(0,0) 



Figure 17 
Cross- 
correlations - 
Longitudinal, 
Window Blanks 
3/4/5, Wide 
Band, Mach 1.6 



Figure 18 
Cross- 
correlations - 
Longitudinal, 
Window Blanks 
3/4/5, Wide 
Band, Mach 
1.95 
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Figure 19 Convection Velocity as a Function of Separation Distance and Mach Number, 

Window Blanks 3/4/5 

5. 2. 1.3 Cross-Correlations at Rear Fuselage 

Figure 20 and Figure 21 show cross-correlations, with longitudinal separation distance as a pa- 
rameter, for transducers on window blank 7. 

Figure 22 and Figure 23 show cross-correlations, with lateral separation distance as a parameter, 
for transducers on window blank 7. 
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Figure 20 
Cross- 
correlations - 
Longitudinal, 
Window 
Blank 7, Wide 
Band, Mach 
0.74 


Figure 21 
Cross- 
correlations - 
Longitudinal, 
Window 
Blank 7, Wide 
Band, Mach 
1.95 
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Figure 22 
Cross- 
correlations - 
Lateral, Window 
Blank 7, Wide 
Band, Mach 
0.74 


Figure 23 
Cross- 
correlations - 
Lateral, Window 
Blank 7, Wide 
Band, Mach 
1.95 
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5.2.2 Cross-Spectra 

The cross-spectral density is defined as the Fourier transform of the cross-correlation. In Section 
5.2.1, the cross-correlation is defined with normalization. Here, we define the cross-spectrum with- 
out normalization: 


oo \ T 

®( x ,^/) = J p(*,t)p(x+z=,t+T)dt 

—CD T— >00 _rp 


z lmiT dx 


where /is frequency. Again, the location vector x is omitted below and is implied by context (win- 
dow blank). 

As shown in Figure 2, window blanks 1 and 7 only offer three longitudinal and two lateral separa- 
tion distances. The cluster of transducers in and between window blanks 3 to 5 offer a much 
greater range of longitudinal separations, but only two in the lateral direction. 

All cross-spectra have dimensions Pascal-squared per Hertz. Ref. 13 was used to compute the 
predicted values. 


5.2.2 . 1 Cross-Spectra at Forward Window Blank 

Figure 24 through Figure 29 contain measured and predicted cross-spectral values as a function of 
longitudinal separations with frequency as the parameter, one set for a mach number of 0.74, an- 
other set for mach 1.95. Measured magnitudes are generally lower than predicted ones. At the 
subsonic Mach number, decay rates with distance are comparable. At the supersonic speed, the 
measured cross spectra show little decay. Phases compare well in all cases. 

Figure 30 and Figure 31 show cross-spectral comparisons (magnitude only) for laterally separated 
transducers. Again, the measured magnitudes and decay rates are lower than the predicted ones. 
At Mach 1.95, predicted cross-spectral values coincide for all frequencies; the measured values 
tend to cluster closer together than at Mach 0.74. 
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Figure 24 Cross-Spectra - Longitudinal, Window Blank 1, Measured and Predicted at Mach 0.74, 32 to 100 Hz 
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Phase(<j)), Radians 


Measured ())(c. 1 ,0,f=const) M=0.74 Alt=1 6400 ft Predicted ,0,f=const) M=0.74 Alt=1 6400 ft 






Figure 25 Cross-Spectra - Longitudinal, Window Blank 1, Measured and Predicted at Mach 0.74, 125 to 500 Hz 
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Phased), Radians 


Measured 4>(£, 1 ,0,f=const) M=0.74 Alt=16400 ft 



Predicted 4>(£, 1 ,0,f=const) M=0.74 Alt=1 6400 ft 





Figure 26 Cross-Spectra - Longitudinal, Window Blank 1, Measured and Predicted at Mach 0.74, 600 to 2500 Hz 
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Figure 27 Cross-Spectra - Longitudinal, Window Blank 1, Measured and Predicted at Mach 1.95, 32 to 100 Hz 
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Phased), Radians 


Measured ,0,f=const) M=1.95 Alt=55400 ft 


Predicted ,0,f=const) M=1.95 Alt=55400 ft 





Figure 28 Cross-Spectra - Longitudinal, Window Blank 1, Measured and Predicted at Mach 1.95, 125 to 500 Hz 
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Phased), Radians 


Measured ,0,f=const) M=1 .95 Alt=55400 ft 



Predicted <|>(!; ,0,f=const) M=1 .95 Alt=55400 ft 





Figure 29 Cross-Spectra - Longitudinal, Window Blank 1, Measured and Predicted at Mach 1.95, 600 to 2500 Hz 
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Figure 30 Cross-Spectra - Lateral, Window Blank 1, Measured and Predicted at Mach 0.74 
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Figure 31 Cross-Spectra - Lateral, Window Blank 1, Measured and Predicted at Mach 1.95 
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5. 2. 2. 2 Cross-spectra at Center Window Blank Cluster 

Figure 32 through Figure 46 show cross-spectral values plotted as a function of longitudinal 
(stream-wise) separation distance obtained at window blanks 3 through 5, and at transducers SI 
and S2 between windows. The measured narrow band cross spectral real and imaginary parts 
were averaged over third octave bands. Measured and predicted (using Efimtsov’s model) magni- 
tude and phase are shown for three different frequency regions on three separate graphs for each 
of five Mach numbers. We chose not to normalize the data (by dividing by the rms value for the 
selected frequency band) in order to provide magnitude comparisons between frequencies and be- 
tween Mach numbers. Here are some observations: 

• Comparisons of measured with predicted data is more favorable at Mach numbers away from 
Mach 1 (0.74 and 1.95) than for Mach 0.88 and 1.2. 

• Predicted magnitude is usually higher than measured at smaller separations. Regarding decay 
rates, predicted magnitude decreases with increasing separation faster than measured. 

• The Mach 1.6 case (starting at Figure 41) shows abnormally high cross-spectral values. The 
measured low frequency phase behaves much more like the predicted than at other Mach 
numbers. We have no explanation for this occurrence. 

• Measured cross spectral values often peak around 500 Hz for all separation distances. This is 
in agreement with predictions. 

• The measured phase behaves much more erratically than the predicted one in the low fre- 
quencies. In the mid frequencies, measured and predicted phases often compare well. As one 
would expect, phase behavior at the higher frequencies is hard to discern. 

Figure 47 and Figure 48 show cross-spectral values (magnitude only) plotted as a function of lat- 
eral (cross-stream) separation distance obtained at window blank 4. Only two such separations 
were available: 60 and 120 mm. Some observations: 

• At Mach 0.74, the magnitude is much smaller than predicted, whereas it is comparable at Mach 
2 . 

• The slopes with separation distance are generally comparable to the prediction for Mach 0.74, 
but not for Mach 1.95. 

• At Mach 1.95, the prediction almost independent of frequency; the measured data shows a 
definite frequency dependence. 

• Measured and predicted data both exhibit a trend towards lower values with increasing Mach 
number. 

No attempt was made to improve upon the cross functional model. 
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Figure 32 Cross-Spectra - Longitudinal, Central Window Blanks, Measured and Predicted at Mach 0.74, 32 to 100 Hz 
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Phase(4>), Radians 


Measured ,0,f=const) M=0.74 Alt=1 6400 Predicted <|)(^ 1 ,0,f=const) M=0.74 Alt=1 6400 






Figure 33 Cross-Spectra - Longitudinal, Central Window Blanks, Measured and Predicted at Mach 0.74, 125 to 500 Hz 
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Phased), Radians 


Measured ,0,f=const) M=0.74 Alt=16400 



Predicted 0(^,0, f=const) M=0.74 Alt=16400 





Figure 34 Cross-Spectra - Longitudinal, Central Window Blanks, Measured and Predicted at Mach 0.74, 600 to 2500 Hz 
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Figure 35 Cross-Spectra - Longitudinal, Central Window Blanks, Measured and Predicted at Mach 0.88, 32 to 100 Hz 
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Figure 36 Cross-Spectra - Longitudinal, Central Window Blanks, Measured and Predicted at Mach 0.88, 125 to 500 Hz 
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Figure 37 Cross-Spectra - Longitudinal, Central Window Blanks, Measured and Predicted at Mach 0.88, 600 to 2500 Hz 
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Figure 38 Cross-Spectra - Longitudinal, Central Window Blanks, Measured and Predicted at Mach 1.2, 32 to 100 Hz 
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Figure 39 Cross-Spectra - Longitudinal, Central Window Blanks, Measured and Predicted at Mach 1.2, 125 to 500 Hz 
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Measured (|)(^ 1 ,0,f=const) M=1.2 Alt=41000 



Predicted 0(^,0, f=const) M=1.2 Alt=41000 





Figure 40 Cross-Spectra - Longitudinal, Central Window Blanks, Measured and Predicted at Mach 1.2, 600 to 2500 Hz 
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Figure 41 Cross-Spectra - Longitudinal, Central Window Blanks, Measured and Predicted at Mach 1.6, 32 to 100 Hz 


44 


Phased), Radians 


Measured 0(^,0, f=const) M=1.6 Alt=45000 



^ , mm 


Predicted (|>(5 1 ,0,f=const) M=1.6 Alt=45000 



^ , mm 




^ , mm , mm 

Figure 42 Cross-Spectra - Longitudinal, Central Window Blanks, Measured and Predicted at Mach 1.6, 125 to 500 Hz 
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Figure 43 Cross-Spectra - Longitudinal, Central Window Blanks, Measured and Predicted at Mach 1.6, 600 to 2500 Hz 
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Figure 44 Cross-Spectra - Longitudinal, Central Window Blanks, Measured and Predicted at Mach 1.95, 32 to 100 Hz 
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Figure 45 Cross-Spectra - Longitudinal, Central Window Blanks, Measured and Predicted at Mach 1.95, 125 to 500 Hz 
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Figure 46 Cross-Spectra - Longitudinal, Central Window Blanks, Measured and Predicted at Mach 1.95, 600 to 2500 Hz 
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Figure 47 Cross-Spectra - Lateral, Central Window Blanks, Measured and Predicted at Mach 0.74 
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Figure 48 Cross-Spectra - Lateral, Central Window Blanks, Measured and Predicted at Mach 1.95 
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5. 2. 2. 3 Cross-Spectra at Aft Window Blank 

Figure 49 through Figure 54 contain cross-spectral values for longitudinal separations. At Mach 
0.74, measured and predicted values (magnitude and phase) are comparable; measured decay 
rates are shallower. At Mach 1 .95, the measured low frequency magnitudes for the smallest sepa- 
ration (40mm) are exceptionally high. At higher frequencies the range of values is much wider 
than the predicted ones. Phases compare very well except at the low frequencies. 

Figure 55 and Figure 56 show cross-spectral values (magnitude only) for lateral separations. At 
Mach 0.74, the low frequency magnitudes and decay rates are comparable; high frequency meas- 
ured magnitudes are significantly lower. As noted before, at Mach 1.95, the prediction is almost 
independent of frequency in contrast to the measurements. 
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Figure 49 Cross-Spectra - Longitudinal, Window Blank 7, Measured and Predicted at Mach 0.74, 32 to 100 Hz 
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Figure 50 Cross-Spectra - Longitudinal, Window Blank 7, Measured and Predicted at Mach 0.74, 125 to 500 Hz 
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Figure 51 Cross-Spectra - Longitudinal, Window Blank 7, Measured and Predicted at Mach 0.74, 600 to 2500 Hz 
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Figure 52 Cross-Spectra - Longitudinal, Window Blank 7, Measured and Predicted at Mach 1.95, 32 to 100 Hz 
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Figure 53 Cross-Spectra - Longitudinal, Window Blank 7, Measured and Predicted at Mach 1.95, 125 to 500 Hz 
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Figure 54 Cross-Spectra - Longitudinal, Window Blank 7, Measured and Predicted at Mach 1.95, 600 to 2500 Hz 
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Figure 55 Cross-Spectra - Lateral, Window Blank 7, Measured and Predicted at Mach 0.74 
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Figure 56 Cross-Spectra - Lateral, Window Blank 7, Measured and Predicted at Mach 1.95 
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6. Efimtsov Auto Spectrum Model Adjustments 

Comparisons of the measured flight data and the predictions from the Efimtsov model showed two 
characteristics which the Efimtsov model did not predict. The first was a broad band spectral peak 

2 k f 5* 

around a Strouhal number of 0.6 where the Strouhal number is calculated as — ; . The meas- 

U 0 

ured data also showed a slightly steeper roll-off at high frequencies (above 1 kHz) than the pre- 
dicted roll-off. 

6.1 Broad Peak 

The broad band peak’s magnitude is of the order of 2 dB. In some of the data it is obvious and 
prominent, in some of the data it is non-existent. Here are arguments for the presence of the 
broad peak: 

• Intuitively, one would expect that certain frequency regions would contribute more strongly to 
turbulence energies according to the length scales imposed on the flow by the boundary layer 
thickness; such behavior has been shown to exist for free shear layers in jet plumes. The 
measured data’s behavior supports this reasoning by exhibiting its spectral peak at a constant 
Strouhal number. 

• Schewe [Ref. 16] shows a figure (here reproduced as Figure 57) with spectra that exhibit the 
broad spectral peak. 



Figure 57 Excerpt from Ref. 16: ‘Corrected power spectral densities 

• The scatter in the low frequency auto spectral data can be attributed - at least in part - to trans- 
ducers not being completely flush mounted. Efimtsov’s [Ref. 8] wind tunnel data indicates that 
the probability is higher of measuring a higher level due to flushness deviations than one below 
the ‘true’ value. One might therefore believe the lower portion of the scattered points more, 
leading to the broad peak asserting its presence. 
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• The low frequency data for window blanks in the back of the aircraft may be contaminated by 
engine noise, air as well as structure borne, leading to higher measured levels than would be 
measured in a ‘pure’ boundary layer. 

Here are arguments against the presence of the broad peak: 

• The broad peak’s magnitude is small; it is within the scatter of the data. 

• The majority of past models exhibit no broad peaks (the exception is the ESDU model). 

To the authors, the arguments ‘for’ the peak were more compelling. Therefore, a procedure was 
developed to add an appropriate adjustment to Efimtsov’s model. The peak occurs in a frequency 
region which is weighted strongly in human response to noise metrics for interior noise; it is there- 
fore important to model loads there as accurately as possible. 

6.2 High Frequency Roll-Off 

The Efimtsov model predicts a much shallower high frequency negative slope than all other mod- 
els, as well as the TU-144LL data. We decided to adjust the Efimtsov model to follow the TU- 
144LL after correcting for transducer size (see Section 4, ‘Data Correction’). 

6.3 Efimtsov Adjustment Functions 

Two analytical functions were used to adjust the predicted curves to more accurately reproduce the 
measured SPL spectra. 

The broadband peak was modeled by a log-normal distribution equation: 


Ciexp(-(log(*)-log(* 1 ))"j 


where St is the Strouhal number, and ST is a reference value. 


To account for the slope of the high frequency roll-off, another function is added to the prediction 
equation. The Efimtsov model over-predicts for flight conditions below Mach 1.65 and under - 
predicts for conditions above 1 .65 which leads to the following adjustment function: 

(M -1.65)log 10 (/) 

Since only the high frequency slope needs adjustment, a tanh function centered at 1000 Hz is in- 
cluded. The final adjustment equation is: 


®(/) 


adjusted 


-9(f) 


predicted 


2.5 exp 


f r 

In 


2 k8‘ 


V V u. 


-(/)]- *n(0.6) 

J 


2 \ 


+ 


tanh 


log 


10 ' '^1000 


+ 1 (M -1.65)log 10 (/) 


Figure 58 shows an example of this function. 
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Frequency, Hz 


Figure 58 Sample Adjustments for Mach 2, S* = \ ft , U 0 =1800 fps. 

Figure 59 shows the original Efimtsov model and the adjusted model for Mach 1.95 55,400ft flight 

condition as well as measured flight test data. Other conditions are shown in Figure 60 - Figure 

64. 

6.4 Comparisons of Adjusted Efimtsov Prediction to Measured TU-144LL Data 

Some observations: 

• The adjusted model follows the measured data better on forward and center window blanks 
than on the aft. 

• The high frequency roll-off of the data is well modeled, except that the model does not follow a 
data trend where further aft sensors exhibit a greater high frequency slope, crossing over the 
spectra of further forward sensors. 

• At low frequencies, no adjustments were made to the Efimtsov model because there was no 
trend to be found in the widely scattered data. 
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Predicted and Measured SPL Data, Mach 1.95, Alt=55,400 ft 
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Measured WB 7 
Efimtsov Predicted WB 7 
Efimtsov Adjusted WB 7 
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Efimtsov Predicted WB 4 
Efimtsov Adjusted WB 4 
Measured Forward 
Efimtsov Predicted Forward 
Efimtsov Adjusted Forward 
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Frequency, Hz 

Figure 59 Adjusted Efimtsov Model and In-Flight SPL Data (M 1.95, 55,400 ft, Window 

Blanks 1, 4, and 7) 
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Predicted and Measured SPL Data, Mach 0.74, Alt=1 6400 ft 



Measured Window Blank 7 
Predicted Window Blank 7 
Measured Window Blank 4 
Predicted Window Blank 4 
Measured Window Blank 1 
Predicted Window Blank 1 
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Frequency, Hz 


Figure 60 Adjusted Efimtsov Model and Tu144 In-Flight SPL Data (M 0.74, 16,400 ft, Win- 
dow Blanks 1, 4, and 7) 
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O(f), dB, (ref: 2x10 Pa) 


Predicted and Measured SPL Data, Mach 0.877, Alt=29500 ft 



Measured Window Blank 7 
Predicted Window Blank 7 
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Predicted Window Blank 4 
Measured Window Blank 1 
Predicted Window Blank 1 
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Figure 61 Adjusted Efimtsov Model and Tu144 In-Flight SPL Data (M 0.88, 29,500 ft, Win- 
dow Blanks 1, 4, and 7) 
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Predicted and Measured SPL Data, Mach 1 .2, Alt=41000 ft 



Figure 62 Adjusted Efimtsov Model and Tu144 In-Flight SPL Data (M 1.2, 41,000 ft, Window 

Blanks 1, 4, and 7) 
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O(f), dB, (ref: 2x10 Pa) 


Predicted and Measured SPL Data, Mach 1.6, Alt=45000 ft 
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Figure 63 Adjusted Efimtsov Model and Tu144 In-Flight SPL Data (M 1.6, 45,000 ft, Window 

Blanks 1, 4, and 7) 
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<t(f), dB, (ref: 2x10“°Pa) 


Predicted and Measured SPL Data, Mach 1 .95, Alt=55400 ft 



Figure 64 Adjusted Efimtsov Model and Tu144 In-Flight SPL Data (M 1.95, 55,400 ft, Win- 
dow Blanks 1, 4, and 7) 
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6.5 Comparison to Other Flight Test Data 

A comparison of other measured data could be used in order to determine if the adjustments are 
specific to the NASA-Boeing Tu144LL test or if the adjustments could be applied to all flight vehi- 
cles. The existing supersonic flight data is limited to the Concorde [Ref. 9], and XB-70 [Ref. 8], 
Several sources of wind tunnel data exist but will likely be of limited use for direct comparison pur- 
poses. 

Figure 65 shows the adjusted Efimtsov predicted SPL and those measured on the Concorde [Ref. 
9], A broad spectral peak seems to be present in the Concorde data as well. The gross outlier in 
window blank 2 probably needs to be attributed to a flushness problem. If the Corcos correction 
were to be applied to the Concorde data the measured and predicted high frequency roll-offs would 
match a lot better. 
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Concorde Predicted and Measured SPL Data, Mach 1.98, Alt=52300 ft 
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Predicted Window Blank 1 


Figure 65 Adjusted Efimtsov Predicted SPL and Measured Flight Test Data from Concorde 

(Mach 2) 


71 


7. Consequences of Model Adjustments to HSCT Interior Noise 

What follows is a qualitative assessment of what changes one might expect in HSCT interior noise 
predictions that are based on the Efimtsov model: 

• The broad spectral peak is added in a frequency range roughly from 200 to 3000 Hz. In the 
higher end of this range the human ear is particularly sensitive. In the lower end of that range, 
fuselage structures often exhibit coincidence effects (particularly honeycomb) or panel reso- 
nances (particularly skin/stringer). 

• The small change of the high frequency roll-off slope decreases the high frequency portion at 
Mach numbers below 1.65, and increases it above that Mach number. 

Without the benefit of a detailed quantitative analysis, one might engage in educated speculation 
that high Mach number cruise interior levels will increase somewhere between 1 to 3 dB. 

For most of the length of the fuselage (middle and aft) the model adjustment adds noise energy 
mostly at frequencies below 1000 Hz, regardless of Mach number. This is the area where active 
noise control shows the most promise. 

8. Recommendations for Future Work 

The following items should be considered for future work: 

• Auto-Spectra: 

• Use Boeing data (777, MD-90) to further evaluate subsonic performance of the adjusted 
model. 

• Explanation of higher than currently predicted levels in the aft portion of the fuselage, 
found in Concorde as well as TU-144 measurements. Possibly further adjust the 
model. 

• Use A. Frendi’s numerical prediction method to investigate the high Reynolds number 
and the thick boundary layer situation on the aft fuselage. 

• Include in the model effects of pressure gradients on the boundary layer, particularly the 
negative gradient on the tail cone section. 

• Cross-Spectra: 

• Correct the measured cross-functions for the influence of local flow direction (measured 
on the TU-1 44 in 1 999 [Ref. 1 8]) 

• Adjust the cross-functional model using measured cross-functions using more of the 
TU-144 data than presented here 

• Use ideas from Smol’yakov-Tkachenko and further develop a model for high Reynolds and 
Mach numbers without the assumption of multiplicative separability of the functions of the time 
and space separation variables. 
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10. Appendix A - ESDU Equations 

Equations for calculating fluctuating pressure levels derived by ESDU [Ref. 3]: 

L a =178.441947 + log 10 (M)(39.29984 + log 10 (M)(3.0984996-log 10 (M)(25.32020 + 25. 18335 -log 10 (M)))) 
L b = -5.2735 -log 10 (M )(0. 19002 + log 10 (M)(2.29509-log 10 (M)(4.83039 + 5.93327 -log 10 (M )))) 

=0. 16237-log 10 (M)(0.00510-log 10 (M)(0.01416-log 10 (M)(0.30915 + 0.37894 -log 10 (M)))) 


L = 20 x log. 


1481.34 

V n ^ 


A = A, + log 10 (Re) X (4 + log 10 (Re) *L C )-L p 


4=10xlog 10 -f 


A =-8.24-logio(*) x 

(l2.84 + log 10 (*)(6.95 + log 10 (*)(24.75 + log 10 (*)(29.30 + log 10 (*)(l4.26 + log 10 (*)(3.17 + 0.27-log 10 (*)))) 


® (f t ) = A + A “■ A, + 10 x lo §io (bandwidth (f , )) 
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11. Appendix B - Fourier Transform of a Turbulence Cross-Correlation 
Model 


This work was performed with the expectation of creating an improved cross-spectral pressure 
fluctuation model, but could not be completed due to lack of time. The partial effort is documented 
here for the time when it can be continued. 


When computing cross functions using the TU-144 data, cross-correlations turn out to be nicely 
smooth curves to which a model can be fitted with little uncertainty. Cross-spectra are harder to fit 
since they contain more noise and consist of real and imaginary parts. The approach to modeling 
cross-spectra is therefore to model the cross-correlations and to Fourier transform that model to 
obtain the cross-spectrum. 

One way of modeling a cross-correlation function as a function of time delay r is shown in Eq. 1. 
This form was chosen such as to be able to fit measured correlations with functions that have rela- 
tively simple Fourier transforms. Ref. 12 provides Fourier transform formulas. 
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Eq. 1 


Meaning of Symbols: 

r = time delay of peak of cross-correlation 
z si = time delay scale factors 
r = magnitude constants 
f 0 = frequency of spectral peak 

The objective is to find the Fourier transform of this function, defined by 
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Eq. 3 


75 



Each element on the right hand side of Eq. 3 can be Fourier transformed separately: 

F { r o} =r o^(f - 0 ) 


Eq. 4 
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